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(N ' We review recent results on the inclusive electroproduction of light hadrons at 

ncxt-to-leading order in the parton model of quantum chromo dynamics impie- 
ty ~ ) , mented with fragmentation functions and present updated predictions for HERA 
' experiments based on the new AKK set. We also discuss phenomenological impli- 
es ~ ) ' cations of these results. 

■ , 1. Introduction 

Oh. 

In the framework of the parton model of quantum chromodynamics (QCD), 
the inclusive production of single hadrons is described by means of frag- 
mentation functions (FFs) D%(x,n). At lowest order (LO), the value of 
. D%(x, ix) corresponds to the probability for the parton a produced at short 

| distance to form a jet that includes the hadron h carrying the frac- 

tion x of the longitudinal momentum of a. Analogously, incoming hadrons 
and resolved photons are represented by (non-perturbative) parton density 
functions (PDFs) F^(x,^). Unfortunately, it is not yet possible to calcu- 
late the FFs from first principles, in particular for hadrons with masses 
smaller than or comparable to the asymptotic scale parameter A. How- 
ever, given their x dependence at some energy scale /i, the evolution with 
fi may be computed perturbatively in QCD using the timclikc Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. Moreover, the factor- 
ization theorem guarantees that the D^(x,ijl) functions are independent of 
the process in which they have been determined and represent a universal 
property of h. This entitles us to transfer information on how a hadronizes 
to h in a well-defined quantitative way from e + e~ annihilation, where the 
measurements are usually most precise, to other kinds of experiments, such 
as photo-, lepto-, and hadroproduction. Recently, FFs for light charged 
hadrons with complete quark flavour separation were determined 1 through 
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Figure 1. Parton-model representation of ep — > eh + X. 

a global fit to e + e~ data from LEP, PEP, and SLC including for the first 
time the light-quark tagging probabilities measured by the OPAL Collabo- 
ration at LEP, 2 thereby improving previous analyses. 3,4 

The QCD-improved parton model should be particularly well applica- 
ble to the inclusive production of light hadrons carrying large transverse 
momenta (jpt) m deep-inelastic lepton-hadron scattering (DIS) with large 
photon virtuality (Q 2 ) due to the presence of two hard mass scales, with 
Q 2 ,Pt ^> A 2 . In Fig. 1, this process is represented in the parton-model pic- 
ture. The hard-scattering (HS) cross sections, which include colored quarks 
and/or gluons in the initial and final states, are computed in perturba- 
tive QCD. They were evaluated at LO more than 25 years ago. 5 Recently, 
the next-to-leading-order (NLO) analysis was performed independently by 
three groups. 6,7,8 A comparison between Refs. 7, 8 using identical input 
yielded agreement within the numerical accuracy. 
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The cross section of e + p — > e + 7r° + X in DIS was measured in various 
distributions with high precision by the HI Collaboration at HERA in the 
forward region, close to the proton remnant. 9,10 This measurement reaches 
down to rather low values of Bjorken's variable xb = Q 2 /(2P-q), where P 
and q are the proton and virtual-photon four-momenta, respectively, and 
Q 2 = —q 2 , so that the validity of the DGLAP evolution might be challenged 
by Balitsky-Fadin-Kuraev-Lipatov (BFKL) dynamics. 

In Ref. 7, the HI data 9 ' 10 were compared with NLO predictions evalu- 
ated with the KKP FFs. 3 In Sec. 2, we summarize the analytical calculation 
performed in Ref. 7. In Sec. 3, we present an update of this comparison 
based on the new AKK FFs. 1 Our conclusions are summarized in Sec. 4. 

2. Analytical calculation 

The partonic subprocesses contributing at LO are 

7* + q -» q + a, 
7* + « -> g + q, 

J*+g^q + q, (1) 

where q represents any of the nj active quarks or antiquarks and it is 
understood that the first of the final-state partons is the one that fragments 
into the hadron h. 

At NLO, processes (1) receive virtual corrections, and real corrections 
arise through the partonic subprocesses 
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-> q' + q' + q, 



where q' ^ q,q. The virtual corrections contain infrared (IR) singulari- 
ties, both of the soft and/or collincar types, and ultraviolet (UV) ones, 
which are all regularized using dimensional rcgularization with D = 4 — 2e 
space-time dimensions yielding poles in e in the physical limit D — ► 4. The 
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latter arise from one-loop diagrams and are removed by renormalizing the 
strong-coupling constant and the wave functions of the external partons in 
the respective tree-level diagrams, while the former partly cancel in combi- 
nation with the real corrections. The residual IR singularities are absorbed 
into redefinitions of the PDFs and FFs. We extract the IR singularities in 
the real corrections by performing the phase space integrations using the 
dipole subtraction formalism. 11 

3. Comparison with HI data 

We work in the modified minimal-subtraction (MS) renormalization and 
factorization scheme with rif — 5 massless quark flavors and identify the 
renormalization and factorization scales by choosing /i 2 = £[Q 2 + (p^) 2 ]/2, 
where the asterisk labels quantities in the j*p center-of-mass (cm.) frame 
and £ is varied between 1/2 and 2 about the default value 1 to esti- 
mate the theoretical uncertainty. At NLO (LO), we employ set CTEQ6M 
(CTEQ6L1) of proton PDFs, 12 the NLO (LO) set of AKK FFs, 1 and the 
two- loop (one- loop) formula for the strong-coupling constant a{ nf \p) with 
A< 5 ) = 226 MeV (165 MeV). 12 

The HI data 9,10 were taken in DIS of positrons with energy E e = 
27.6 GeV on protons with energy E p = 820 GcV in the laboratory frame, 
yielding a cm. energy of y/S = 2y/E e E p = 301 GeV. The DIS phase 
space was restricted to 0.1 < y < 0.6 and 2 < Q 2 < 70 GeV 2 , where 
y = Q 2 /(xbS). The ir° mesons were detected within the acceptance 
cuts pj, > 2.5 GcV (except where otherwise stated), 5° < 9 < 25°, and 
xe > 0.01, where 9 is their angle with respect to the proton flight direction 
and E — xeE p is their energy in the laboratory frame. The comparisons 
with our updated LO and NLO predictions are displayed in Figs. 2(a)-(d). 
The QCD correction (K) factors, i.e. the NLO to LO cross section ratios, 
are presented in the downmost frame of each figure. 

Comparison of Figs. 2(a)-(d) with Figs. 3, 5(a), 6(c), and 7 in Ref. 7, 
where the KKP FFs 3 were used, reveals that the update of our FFs, from 
set KKP to set AKK, 1 has hardly any visible impact on the theoretical 
predictions considered here. This may be understood by observing that the 
OPAL light-quark tagging probabilities for charged pions, 2 included in the 
AKK analysis, agree well with the assumption made in the KKP one that 
(x, jUo) = (x, /j,q) at the starting scale p, of the DGLAP evolution. 
In Figs. 3(a) and (b), 8 the HI data 10 on da/dp* T for 2 < Q 2 < 4.5 GeV 2 , 
4.5 < Q 2 < 15 GcV 2 , or 15 < Q 2 < 70 GeV 2 and on da/dx B for p* T > 
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3.5 GeV and 2 < Q 2 < 8 GeV 2 , 8 < Q 2 < 20 GcV 2 , or 20 < Q 2 < 70 GeV 2 , 
respectively, are compared with the LO and NLO predictions evaluated with 
the KKP FFs 3 or those by Kretzer (K). 4 While the LO predictions based 
on the KKP and K sets agree very well, the NLO predictions based on the 
K set appreciably undershoot those based on the KKP set. If it were not 
for the theoretical uncertainty, one might conclude that the HI data prefer 
the KKP set at NLO. 

From the downmost frames in Figs. 2(a)-(d), we observe that the K 
factors are rather sizeable, although the /i values arc reasonably large. In 
Fig. 4, 8 the impact of the HI forward-selection cuts on the K factor is 
studied for the case of da/dx B for 2 < Q 2 < 8 GeV 2 and p* T > 3.5 GeV. 
Towards the lower end of the considered xb range, the K factor reaches one 
order of magnitude if these cuts are imposed [see also Fig. 2(c)]. However, 
if the latter are removed, the K factor collapses to acceptable values of 
around 3. From this finding, we conclude that these cuts almost quench 
the LO cross section. In other words, in the extreme forward regime, the 
latter is effectively generated by the 2^3 partonic subprocesses of Eq. (2). 

It is interesting to investigate the relative importance of the tagged par- 
tons, i.e. the one (a) that originates from the proton and the one (b) that 
fragments into the hadron. In Fig. 5, 8 the NLO contributions from the 
four most important ab channels to da/dxs for 2 < Q 2 < 8 GeV 2 and 

> 3.5 GeV with the HI forward-selection cuts are shown together with 
the total LO contribution. We observe that the gg channel makes up ap- 
proximately two thirds of the cross section in the 1ow-x_b regime. 

4. Conclusions 

We calculated the cross section of ep — ► e-zr + X in DIS for finite values of 
p* T at LO and NLO in the parton model of QCD 7 using the new AKK FFs 1 
and compared it with a precise measurement by the HI Collaboration at 
HERA. 9 ' 10 

We found that our LO predictions always significantly fell short of the 
HI data and often exhibited deviating shapes. However, the situation dra- 
matically improved as we proceeded to NLO, where our default predictions, 
endowed with theoretical uncertainties estimated by moderate unphysical- 
scale variations, led to a satisfactory description of the HI data in the 
preponderant part of the accessed phase space. In other words, we en- 
countered K factors much in excess of unity, except towards the regime of 
asymptotic freedom characterized by large values of p* T and/or Q 2 . This was 
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unavoidably accompanied by considerable theoretical uncertainties. Both 
features suggest that a reliable interpretation of the HI data within the 
QCD-improved parton model ultimately necessitates a full next-to-next-to- 
leading-order analysis, which is presently out of reach, however. For the 
time being, we conclude that the successful comparison of the HI data 
with our NLO predictions provides a useful test of the universality and the 
scaling violations of the FFs, which are guaranteed by the factorization 
theorem and are ruled by the DGLAP evolution equations, respectively. 

Significant deviations between the HI data and our NLO predictions 
only occurred in certain corners of phase space, namely in the photopro- 
duction limit Q 2 — > 0, where resolved virtual photons are expected to con- 
tribute, and in the limit r\ — > oo of the pseudorapidity r\ = — ln[tan(#/2)], 
where fracture functions are supposed to enter the stage. Both refinements 
were not included in our analysis. Interestingly, distinctive deviations could 
not be observed towards the lowest x b values probed, which indicates that 
the realm of BFKL dynamics has not actually been accessed yet. 
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Figure 2. HI data on (a) da/dp^ and (b) da/dx E for 2 < Q 2 < 4.5 GeV 2 , 4.5 < Q 2 < 
15 GeV 2 , or 15 < Q 2 < 70 GeV 2 , on (c) da/dx B forpj, > 3.5 GeV and 2 < Q 2 < 8 GeV 2 , 
8 < Q 2 < 20 GeV 2 , or 20 < Q 2 < 70 GeV 2 , and on (d) da/dQ 2 from Refs. 9 (open 
circles) and 10 (solid circles) are compared with our default LO (dashed histograms) and 
NLO (solid histograms) predictions including theoretical uncertainties (shaded bands). 
The K factors are also shown. 



February 2, 2008 5:11 



Proceedings Trim Size: 9in x 6in 



kniehl 



700 
600 
500 
400 
300 
200 
100 

LOO 
80 
60 
40 
20 


20 

10 



' 2 < Q 2 < 8 GeV 2 



HI data '04 

LO 

NLO 



. i - 



i < Q < 20 GeV 




H h 



-H- 

20 '< Q" < 70 GeV" 



I 



2 < Q" < 8 GeV" 

8 < Q 2 < 20 GeV 2 

20 < Q 2 < 70 GeV 2 



(d) 




Figure 2. Continued. 
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Figure 3. HI data 10 on (a) da/dp* T for 2 < Q 2 < 4.5 GeV 2 , 4.5 < Q 2 < 15 GeV 2 , 
or 15 < Q 2 < 70 GeV 2 and on (b) da/dx B for p* T > 3.5 GeV and 2 < Q 2 < 8 GeV 2 , 
8 < Q 2 < 20 GeV 2 , or 20 < Q 2 < 70 GeV 2 are compared with the LO and NLO 
predictions evaluated with the KKP 3 or K 4 FFs (taken from Ref. 8). 
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Figure 4. K factors of da/dx B for 2 < Q 2 < 8 GeV 2 and p* T > 3.5 GeV with and 
without the HI forward-selection cuts 10 (taken from Ref. 8). 




Figure 5. Total LO contribution and NLO contributions from the four most important 
ab channels, where a and b are the partons connected with the PDFs and FFs, respec- 
tively, to da/dxg for 2 < Q 2 < 8 GeV 2 and pj, > 3.5 GeV with the HI forward-selection 
cuts 10 (taken from Ref. 8). 



